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0. FOREWORD 

0.1 This Indian Standard (Part 4) was adopted by the Indian Standards 
Institution on 14 June 1985, after the draft finalized by the Reliability 
of Electronic and Electrical Components and Equipment Sectional Com- 
mittee had been approved by the Electronics and Telecommunication 
Division Council. 

0.2 This standard (Part 4) is one of the series of Indian Standards for 
equipment reliability testing and covers procedures for determining point 
estimates and confidence limits from equipment reliability determination 
tests. Other standards available in this series are given on page 28. 

0.3 This standard is largely based on TEC Doc: 56 (Central Office) 92 
Draft, IEC Pub 605 Equipment reliability testing: Part 4 Procedures for 
determining poirt estimates and confidence limits from equipment reli- 
ability determination tests' issued by the International Electrotechnical 
Commission (IEC). 

0.4 This standard is one of a number of Indian Standards on reliability 
of electronic and electrical components and equipment. A list of 
standards published so far in this area is given on page 28. 



1. SCOPE 



1.1 This standard (Part 4) gives recommended numerical and graphical 
methods for determining point estimates and confidence limits of reli- 
ability characteristics from equipment reliability determination tests. 
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2. TERMINOLOGY 

2.0 For the purpose of this standard, the terras and definitions as gives 
in IS : 1885 (Part 39)-1979* and IS: 7690-1975! shall apply. 

3. LIST OF SYMBOLS 

3.1 The following symbols are used in this standard. The 'determina- 
tion point' is the point in time or number of trials at which the point 
estimates or confidence limits are determined. 

a — location parameter in a Weibull distribution; 

b — charateristic life (or scale parameter) in a Weibull 

distribution; 

f(t) — probability density function of time to failure; 

F (t) — cumulative distribution of times to failure, proba- 

bility of failure within time t; 

Fp ( v i> v z) — theoretical value of the F-distribution with v l deg- 

rees of freedom in the numerator and -v t degrees 
of freedom in the denominator at the fractile 
of order p; 

i = order number of a failure based on the time to 

failure (ft); 

k — shape parameter of a Weibull distribution; 

m = true mean time between failures; 

m = point estimate of mean time between failures 

(observed value); 

m p = true mean time to failure, mean life; 

m = point estimate of mean time to failure (observed 

value; 

m (0,t ) = true mean time between failuies for a period of 

time (O, to); 



•Electrotechnicat vocabulary: Part 39 Reliability of electronic and electrical 
items (first revision). 

■^Mathematical guide to the terms and definition; for reliability of electronic 
equipment and components (or parts) used therein. 



IS : 8161 ( Part 4 ) - 1^8$ 

a, t _ _ . 

»i (O, to) ~ point estimate 01 mean time between failures for 

the period of time (0, t ) (observed value); 

n — relevant number of test items or trials; 

^50 (ti) ~ 50 percent (median) ranks of fi (see Table 1); 

r = total number of relevant failures on which the 

determination test is based; 

R = true success ratio; 

Jl (r)=l — F (r) = probability of success; 

f = time (or equivalent, such as distance, cycles or 

other quantites as appropriate); 

J* m = predetermined test time for the mth item; 

t\ = relevant test time of magnitude order number /, 

recorded between the point of time the corres- 
ponding item was put on test and the point of 
time of failure of the item; 

h (v) = theoretical value of the Student t distribution 

with v degrees of freedom at the fractile of 
order p; 

T* = accumulated relevant test time up to the deter- 

mination point when this does not coincide with 
a failure; 

TV — relevant test time accumulated for all items up to 

the determination test point when this coincides 
with a failure; 

A = true constant failure rate; 

A = point estimate of failure rate (observed value); 

a = true standard deviation in normal distribution; 

a = point estimate of the standard deviation; and 

X%( v ) = theoretical value of the X 2 distribution with v 

degrees freedom at the fractile of order p. 
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TABLE 1 MEDIAN RANKS ( P 50 IN PERCENT ) 

( Clauses 3.1, 6.3.1, 7.1. 1, 7-1.2 and 7.2.2.1 ) 
Sample Size ( n ) 



Order 
No (/) 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 


500 


292 


20 6 


15 9 


12 9 


109 


94 


8-3 


7-4 


66 


2 




707 


50-0 


38-5 


313 


26-6 


228 


20-1 


179 


16-2 


3 






79-3 


61-4 


50-0 


42-1 


3<v4 


32-0 


286 


258 


4 








84-0 


686 


578 


500 


440 


39-3 


35-5 


5 










87-0 


73-5 


635 


559 


500 


45-1 


6 












890 


77-1 


679 


60-6 


54-8 


7 














90 5 


79-8 


71-3 


64-4 


8 
















91 7 


82.0 


74- 1 


9 


















92 5 


837 


10 




















933 










Sample Size ( 


n) 










Order 
No. (0 










A 












11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


1 


61 


5-6 


51 


48 


4-5 


4*2 


3-9 


37 


3-5 


34 


2 


14-7 


13 5 


125 


11-7 


10 9 


102 


96 


9-1 


86 


8-2 


3 


235 


216 


20-0 


18 6 


17-4 


163 


154 


14-5 


13-8 


131 


4 


32-3 


29'7 


27-5 


25-6 


239 


22-4 


211 


20 


18-9 


180 


5 


411 


37 8 


350 


32-5 


30 4 


285 


26-9 


25-4 


24- 1 


22-9 


6 


500 


459 


42 5 


395 


369 


34-7 


327 


309 


293 


27-8 


7 


58-8 


54'0 


50 


46 5 


^3-4 


40-8 


38-4 


36.3 


34-4 


32-7 


8 


676 


62 1 


574 


534 


50-0 


46-9 


44-2 


41-8 


396 


37-7 


9 


76-4 


70-2 


649 


60 4 


565 


53 


500 


47-2 


44.8 


426 


10 


85-2 


783 


72-4 


67 4 


630 


59 1 


55-7 


52-7 


500 


47-5 


11 


93-8 


86-4 


79-9 


74-3 


695 


65 2 


615 


581 


55 1 


52-4 


12 




94.3 


87-4 


81-3 


760 


71 4 


672 


636 


603 


57-3 


13 






94-8 


88-2 


82 5 


77-5 


73-0 


69 


655 


62-2 


14 








951 


89 


836 


78-8 


74-5 


706 


672 


15 










954 


89-7 


84-4 


79-9 


75-8 


72-1 


16 












95-7 


90 3 


85-4 


81-0 


77 


17 














960 


90 8 


86-1 


81 9 


18 
















960 


91-3 


86-8 


19 


















96-4 


91-7 


20 




















96 5 



Note — For higher values of n and i use formula: 



4. APPLICABILITY 

4.1 The methods of calculation in this standard may be applied at any 
time or after any number of trials of the reliability testing. The more 
information, that is available when the estimates and limits are deter- 
mined, the higher is the precision obtained. Existing data from earlier 
tests or field observations may be acceptable, provided the data are 
sufficiently complete, well establised and applicable to the situation. The 
standard is not applicable to the pooling of data from different sources, 
for example, original data from tests under different conditions. 
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4.2 Wherever 'time' is used in this standard, this variable may be 
replaced by distance, cycles or other quanties as may be appropriate. 

4.3 If a constant failure rale, a mean time to/between failures, or a 
mean failure rate for a given period of time is to be estimated from the 
reliability determination test, the methods in 6 or 7 should be used. 
These methods are applicable for tests based on time. The method to^ 
be used depends on what assumption can be made regarding the failure 
rate dependence on time and the statistical distribution of time to/ 
between failures. 

4.4 In case of an assumption of a constant failure rate characterising an 
exponential distribution of times to/between failures, the method given 
in 6 is to be used. The test items may be repaired and put back on 
test or not repaired. The graphical method (6.3) is, however, not applic- 
able after the first failure of repaired test items. In case of an assump- 
tion of a non-constant failure rate described by a Weibull distribution 
or a normal (Gaussian) distribution of times to failure, the methods 
given in 7 should be used. These methods are not applicable after the 
first failure of repaired test items. 

4.5 If the success ratio is to be estimated from a reliability determina- 
tion test, the methods in 8 are to be used. These methods are based on 
the number of test items or trials, which are classified as either failed 
or non-failed. The method to be used depends on the number of fail- 
ures. If the number of failures is greater than 9 the normal distribution 
may be used, in other cases the binominal distribution shall be used. 

5. FEATURES OF THE POINT ESTIMATES AND CONFIDENCE 
LIMITS 

5.1 A point estimate is a single numerical value to represent the 
unknown true value of a statistical parameter, such as the failure rate. 
Generally, the point estimate considered here is the 'observed' value 
[ as defined in IS : 1885 ( Part 39 )-1979* and IS : 7690-19751 ]. 

5.2 The confidence limits define a confidence interval around the point 
estimate where the interval includes the true value of the parameter 
being estimated with a certain probability, the confidence level. 

5.3 The confidence interval will be narrower when more information is 
available. The information used is the accumulated relevant test time 
and the number of failures for time based tests, and the number of test 



♦Electrotechnical vocabulary: Part 39 Reliability of electrcnx and electrical 
items ( first revision ). 

tMathematrcal guide to the terms and definitions for reliability of electronic 
equipment and components (or parts) used therein. 
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items or trials and the number of failures for tests based on items or 
trials. 

5.4 The confidence interval may be one-sided or two-sided. In case of 
a one-sided confidence interval, a lower or an upper confidence limit is 
given for the characteristic. In case of a two-sided confidence interval, 
both a lower and an upper limit are given. 

5.5 The preferred confidence level is here chosen to be 90 percent. The 
confidence intervals according to this recommendation therfore include 
the true values of the characteristics with 90 percent probability. If 
other confidence levels are desired, the formulae in this standard may 
be used with appropriate values of the test statistics taken from 
statistical tables. 



6. CONSTANT FAILURE RATE 

<s.O If the true failure rate is assumed to be constant, the numerical and 
graphical methods in 6.1, 6.2 and 6.3 may be used for estimating failure 
rate and mean time between failures (for repaired items) or mean time 
to failure (for non-repaired items). The methods are applicable both 
to repair and non-repair situations. These estimates are based on the 
total number of relevant failures and the accumulated relevant test time 
up to the determination point.' 

6.0.1 The validity of a constant failure rate assumption should prefer- 
ably be tested in accordance with IS : 8161 (Part 6)-1983* before point 
•estimates and confidence limits are calculated. 

For application to mean time between failures, M.T.B.F. = m. 
For application to mean time to ( first ) failure, wf or M.T.T.F.: 
«?f = M.T.T.F. = m. 

For application to failure rate A = ~ 

*<».l Time Terminated Tests 

6.1.0 The accumulated relevant test time up to the determination 
point, T*, may be determined in accordance with Appendix A. The 
•same formulae are used for test with and without, replacement or repair. 
The total number of relevant failure, r, up to the determination point 
is counted. 



•Guide for equipment reliability testing : Part 6 Tests for validity of a constant 
'failure rate assumption. 
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6.1.1 Point Estimates 

The point estimate (observed value) of the failure rate is: 



> - ~ 



The point estimate (observed value) of the mean time to or between 
failures is: 



m 



r 



If no failures are observed up to the termination point, that is, 
r = 0, the following point estimate of the failure rate is recommended: 

a = i/3r* 

Too much reliance should not be placed on estimates based upon 
zero or any small number of failures. 

Care should be taken in planning the tests such that T* will be 
sufficiently large in relation to the expected mean time to failure. 

Note — This recommendation is based upon an investigation by E.L. Welker 
and M. Lipow 'Estimating the exponential failure rate from data with no failure 
events' in Proceedings of 1974 Annual Reliability and Maintainability Sym- 
posium, pages 420-427. 

6.1.2 Confidence Intervals 
6.1.2.1 The confidence limits for the true failure rate, a, with the 
confidence level of 90 percent are shown below {see Table 2 for values 
of the X % distribution) . 

a) One-sided confidence interval, upper limit: 

f X 8 . fl0 (2r+2) . XV 9 o(2 ;- + 2) 

b) Two-sided confidence interval: 
- X 2 .o 5 (2r) < A < - ^ 2 . 95 (2r+2) 



or 



2r ^ 2r 



X 2 q.q-, (2r) X- . 9G ( 2r + 2 ) 

IT* IT* 



If no failures are observed only the one-sided confidence limit with 
an upper limit can be defined. 
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TABLE 2 CHI-SQUARE DISTRIBUTION FRACTILES 

( Clauses 6.1.2.1, 6.1.2.2, 6.2.2.1 and 6.2.2.2 ) 



4605 5-991 

7-779 9-488 

10-65 12-59 

1336 15-51 

1598 1831 

1855 2103 

21-06 23 69 

23-54 26-3 

25-99 28-87 

28-41 31-41 

30-81 33-92 

33-2 36-42 

35-56 38-89 

37-92 41-34 

40-26 43-77 

42-57 46-17 

44-88 48-57 

47*19 50-96 

49-5 5336 

51-81 55-76 

5408 58-11 

6317 67-51 

65-42 69-82 

74-4 79-08 

7663 81-37 

85-53 9053 

8774 92-8 

96-58 101-88 

98-78 104-13 

107-57 H.i-15 

109-76 115-39- 

118-5 124-34 

120-65 12653 

129-25 135-3 

131-4 137-5 

140 146.27 

142-15 148-46 

226-02 233-99 



Degrees of 
Freedom v 


*■ Q'OSV ■> 


2 


■103 


4 


•711 


6 


1 635 


8 


2-733 


10 


394 


12 


5-226 


14 


6-571 


16 


7-962 


18 


9-39 


20 


10-85 


22 


12-34 


24 


13-85 


26 


15-38 


28 


1692 


30 


18-49 


32 


20-C9 


34 


21-7 


36 


23-3 


38 


24-91 


40 


26-51 


* 42 


28-16 


50 


34-76 


52 


3645 


60 


43-19 


62 


44-9 


70 


51-74 


72 


5V47 


80 


6039 


82 


6214 


90 


69 13 


92 


70 89 


100 


77-93 


102 


79-74 


110 


86 96 


112 


88-77 


120 


96 


122 


97-81 


200 


168-28 



1-64 +1-29 +1-64 



Note 1 — Linear interpolation of intermediate values is sufficiently accurate. 
Values for v=2/-+2 for various integer r values are included. 

Note 2 — For higher values, use: 

X p ( v )=[ (z+\/2v-l ) ! ]/2 where z is the corresponding percentage of the standard 
normal distribution given at the foot of each column. 

10 
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6.1.2.2 The confidence limits for the true mean time to or between 
failures m, with a confidence level of 90 percent are shown below (see 
Table 2 for values of the X 2 distribution). 

a) One-sided confidence interval, lower limit: 

2r ^ 27"* 



X*o.9t> (2r + 2) X*o. 90 (2r + 2) 

b) Two-sided confidence interval: 

" ~*V, 5 (^ +TT < W < ™ r o . 05 (^r ) 



or 



2r* 2T* 

-< m< 



X 2 . 95 (2/- +2) ^'"^ X* .o 5 (2r) 
If no failures are observed, only the lower limit can be defined. 

The 90 percent confidence limits as a i function of the number of 
failures, r y are shown in Fig. 1 and Table 3. The limits are here 
expressed by the appropriate multiplier times the relevent point estimate 
•A and m. 

Note — The multipliers may be used in the planning oflife tests by finding out 
the approximate number of failures for a given precision of the estimate. The 
required accumulated test time, T*, for the required number of failures is 
approximately T* = y/>>, where the failure rate X has to be assumed from earlier 
experience. 

6.2 Failure Terminated Tests — The accumulated relevent test time up 
to the determination point, T T , is determined in accordance with 
Appendix A. The same formulae are used for tests with and without 
replacement or repair. 

6.2.1 Point Estimates — The point estimate ( observed value ) of the 
failure rate is: 

>s r 
A ~"' Tr 
The point estimate ( observed value ) of the mean time to or between 
failures is: 

7> 
m = r 

Too much reliance should not be placed on estimates based upon a 
small number of failures. 

Care should be taken in planning the tests such that T* will be 
sufficiently large in relation to the expected mean time to failure. 

11 
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MULTIPLIER FOR FAILURE RATE 
1 
5-0- 



ONE-SIDED UPPER LIMIT 




TWO.-SIDEO INTERVAL 



05 



MULTIPLIER FOR FAILURE RATE 
1 

■>o-o 



TWO-SfDEO INTERVAL 




ONE-SIDED LOWER LIMIT 



0-2 



Fig. l the 90 Percent Confidence limits as a Function of 
Number of Failures For Time— terminated tests 
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TABLE 3 


MULTIPLIERS TO OBTAIN CONFIDENCE LIMITS 








( Clause 6.1.2.2 ) 








Multipliers tor MTTF/MTBF 


Multipliers for Failure Rate 

A. 


R 


M(5 


M( 90 


M(95 


M(5 


M(90 


M(95 




percent ) 


percent ) 


percent ) 


percent ) 


percent ) 


perceni y 


1 


19417 


•257 


•21 


•051 


3-889 


4-743 


2 


5625 


•375 


•317 


•177 


2 662 


3-147 


3 


3669 


•449 


•386 


•272 


2-226 


2-584 


4 


2-927 


•5 


•436 


•341 


1-997 


2-288 


5 


2-538 


•539 


•475 


•394 


1-854 


2103 


6 


2-296 


•569 


•506 


•435 


1-755 


1-974 


7 


2 13 


■594 


■532 


•469 


1 681 


1-878 


8 


2 009 


•615 


•554 


■497 


1 624 


1-804 


9 


1 916 


•633 


•573 


•521 


1-578 


l-745 : 


10 


1-843 


•649 


•589 


■542 


1-54 


1 695 


11 


1782 


•662 


•604 


•56 


1-509 


1-655 


12 


1-732 


•674 


•617 


•577 


1-481 


1-62 


13 


1-69 


•685 


•628 


•591 


1-458 


1-59 


14 


1-654 


•695 


•639 


■604 


1 -437 


1-563 


15 


1 622 


•704 


•649 


•616 


1-419 


F53& 


16 


1592 


•713 


•658 


•627 


1-402 


1-517 


17 


1 566 


•72 


•667 


•638 


1-387 


1-498. 


IS 


1-545 


•727 


•674 


•647 


1-375 


1-482 


19 


1 525 


•733 


•681 


•655 


1-363 


1-467 


20 


1-508 


•739 


•688 


•662 


1 351 


1-452 


25 


1 438 


•764 


•716 


•695 


1-308 


1 396 


30 


1-389 


•782 


•737 


•719 


]-277 


1-356- 


35 


1-352 


■797 


•754 


•739 


1-253 


1-325 


40 


1-324 


•809 


■768 


•754 


1-234 


1-303 


45 


1-301 


■819 


•779 


•768 


1-219 


1 282" 


50 


1-283 


•828 


•79 


•779 


1-206 


1 265 


55 


1-264 


•837 


•8 


•79 


1-194 


1-25 


60 


1-249 


•844 


•808 


•8 


1184 


1-237 



5.2.2 Confidence Intervals 

6.2.2.1 The Confidence limits for the true failure rate with a. 
confidence level of 90 percent are shown below ( See Table 2 for values. 
of the X a distribution ): 

a) One-sided confidence interval, upper limit: 

A<4- X?0 - ! ' 1 2 LL r A < X'Vso (2r ) 
2r 



21 r 



b) Two-sided confidence interval: 



A X g o. OB (2r) A ^X* ., 5 .(2r) 
2r ^ ^ 27> 



or 



X 2 ,, 



(2r) 



27> 



< A< 
13 



X V as ( 2r ) 
llr 
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6.2.2.2 The confidence limits for the true mean time to or between 
failures with a confidence level of 90 percent are shown below, ( see 
Table 2 for values of the X 2 distribution ). 

a) One-sided confidence interval, lower limit: 

2r ^ TTr 

w > m v^ — r^ — \~ orm> 



X-o.,o ( 2r ) " XV.,o ( 2r ) 

b) Two-sided confidence interval: 

2r - m 2r 



< m < 



or 



^ a o.o5 ( 2r ) XVo 5 ( 2r ) 



27> 27> 



X'o.m ( 2/" ) ^ r . 05 ( 2r ) 

f <£.3 Graphical Method 

6.3.0 A method is presented for 'estimating failure rate and mean 
time to first failure or mean time between failures, based on the semi- 
log paper. The tests from which the observations are taken need not 
go on until all test items have failed. 

The semi-log paper should be used only for times to first failure 
and for not less than 4 failures. The method will give point estimates 
and may also give an indication of deviations from a constant failure 
rate. 

6-3.1 A number of items, n, are tested and r failures observed. For 
.each failed item, the relevant test time ti is recorded. 

The times ti are ordered in magnitude: 

t x < r a < t< < tr 

On a semi-log paper ( see Fig. 2 ), the ti values are plotted along 
the linear scale and the reciprocal of one minus the fractional median 
rank, /% ( h )/ 100 along the logarithmic scale. Values of median ranks 
may be found from Table 1. If a constant failure rate assumption holds, 
the plotted points will fit well to a straight line passing through the 
point, f i = 0, ratio = 1. In drawing the line, the central points should 
dominate the determiniration of the slope. This being so, then the 
estimate of the mean time to failure equals the t value on the time axis 
that corresponds with the ratio 2"72 on the vertical axis. The estimate 
of the failure rate is then the inverse of the t value. 

14 
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2-7; 




100 9 300 400 

relevant test -time for each item 

Fig. 2 Example of Use of Semi-Log Paper 
1 

1 - Pz« Oi )/ 100 

1 

m = — 



If the plotted points cannot be fitted by a straight line, the failure 
rate may not be constant. In such a case, other methods as described in. 
7 may apply. 
6.3.2 Example 

6.3.2.1 Figure 2 has been drawn with data collected by putting 10 
items in reliability test which was terminated when 8th failure was 
observed. The observation and calculations are given below: 



Failure 


Test Time (t\) 
and i th Failure 


percent 


1 


Order 


l-Pso( *l) 




( hours ) 




100 


1 


29 


66 


107 


2 


29 


162 


T19 




30'6 


25 8 


1*38 


A 


32'4 


35-5 


1'55 


5 


33'0 


45 1 


F82 


6 


35*3 


54'8 


2"21 


7 


36-1 


64*4 


281 


8 


40-1 


74"! 


3-86 



6.3.2.2 The above data is plotted on semi-log paper in Fig. 2 as 
■described in 6.3,1 and observed value for estimated failure rate at the 



ratio value of 2'72 is equal to ~r-r~, - h 

36 5 



0027 4 failures/hour. 



*Obtained from Table 5. 
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7. NON-CONSTANT FAILURE RATE 

7.0 If the assumption of a constant failure rate is not valid and the dis- 
tribution of times to failure follow a Weibull or a normal distribution, 
the methods given below are applicable. The validity of any of these 
distribution assumptions should preferably be tested in accordance with 
IS : 8161 ( Part 8 )* or using the probability papers described in 7.1 
and 7.2. 

7.0.1 The life history of the test items shall be similar when they are 
put on test. The methods of evaluation are applicable only to the times 
to first failure of each item- The test does not, however, need to be 
continued until all items have failed, except for the method in 7.2.1.1. 

7.0.2 For each failed test item, the relevant test time, t\, is determined. 
The time t\ is defined as the relevant test time between the point of time 
when the item i was put on test, t = 0, and the point of failure, t = tu 
The times t\ are ordered in magnitude 

h < t 2 < l- A < < t, 

Using these times and a Weibull probability paper or a normal probability 
paper it is possible to estimate the parameters of the distribution and at 
the same time obtain an indication of how well the times are fitted to 
the assumed distribution. 

7.1 Weibull Distribution 

7.1.0 The cumulative distribution function for the Weibull distribution 
is: 

F(t) = l-e( t ~ a hc; > a, b > o, k > o 

7.1.1 A graphical method for point estimates is presented below- In 
this method, the Weibuil probability paper is used. This paper is so 
designed that the cumulative distribution function of times to failure 
which follow a Weibull distribution, will become a straight line. See 
Fig. 3 which is an example of such a paper and how it can be used. 

If the time scale on the abscissa is not directly applicable, it can be 
transformed by a scale factor 10°, where c is an interger. 

The observed ti values are plotted at the corresponding median rank 
values ( P 5 o ) ( f i ) on the F ( t ) scale. Median rank values can be found in 
Table I. 



♦Guide for equipment reliability testing: Part 8 Tests for validity of a non- 
constant failure rate assumption ( under consideration ). 
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RELEVANT TEST TIME FOR EACH ITEM (t I -— •*• 

Fig. 3 Example of Use of Weibull Probability^Paper 

7.1.2 The example of Fig. 3 assumes a test sample n = 10 when the 
test is terminated after r = 8 failures. Therefore, there are only 8 
plotted points and in the table of median ranks ( Table 1 ) the column 
for n = 10 is used up to and including / = 8 only. 

7.1.3 If a Weibull distribution with a = holds then the plotted 
points can be fitted will by a straight line. If the points follow a concave 
curve seen from the abscissa, the location parameter a has a positive 
value. If the curve is convex then a is negative. For a curve the 
following procedure may be used: 

Step 1: Plot curved data as described above and illustrated in Fig, 3. 

Step 2: Draw two parallel horizontal lines through the extreme 
failure points, and a third horizontal parallel line at the 
linear mid-point- 

Step 3: Draw three lines perpendicular to the ?-axis from the points 
of intersection of the three horizontal lines and the curve. 

Step 4: Identify the lowest, middle, and highest t value as ti, tm 
and fb respectively. 
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Step 5: Use the following formula to calculate: 
^ ( tn — tm) ( tm — ti) 

( /h — .tm ) ( tm — ti) 

A A. 

Step 6: Subtract a ( including sigh of a ) from each h value. 

Step 7: Replot the data; the points should now approximate a 
straight line, a shifted Weibull distribution of failures. 

7.1.3.1 The data plotted in Fig. 3 was collected on 10 items kept on 
reliability test and the test was terminated at 8th failure. The observed 
failures and calculations as described above are given below: 

Failure Test Time, at i th Failure Medium Rank Value 

Order ( hours) percent* 

1 400 6'6 

2 510 162 

3 620 258 

4 720 35'5 

5 820 45*1 

6 920 54-8 

7 1 050 64-4 

8 1 300 74T 

7.1.3.2 From above data, h = 400, t h = 1 300 and tm = 620 and 
hence a = 294"8 ( & 300 ). 

7.1.4 To estimate the shape parameter, k, a line is drawn through the 
centre of the circular scale parallel to the fitted straight line. The value 
of k is then read off the circular scale at the intersection. 

The estimate of the characteristic life or scale parameter, b, is equal 
to the t value of the intersection of the fitted straight line and dashed 
line on the vertical axis at 6396. 

7.1.5 To estimate mean life wf refer to Fig. 4. This gives for different 
k values the value of F( t ) when t = rnF. Then by reference to the 
ordinate scale of F ( t) in Fig. 3, m? can be read off from the original 
plot ( not the straight line derived from a ) at the abscissa scale. 
The original plot may also be used to estimate the fraction F ( t), of 
the population having times to failure less than a given t value. 

7.1.6 Other available types of graph paper, while basically similar, 
may use different methods to estimate the parameters as instructed. 



♦Taken from Table 5. 
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7.2 Normal Distribution 



7.2.0 The probability density function of the normal distribution is: 

This distribution is useful for times to failure, provided that wf > 3a. 

7.2.1 Numerical Method 

1. 1.1.1- Pnirtt 0Vtiwinft>V TT-t£» nilmonVnl n.r*4-',-r^n*-nr, ^f+U„ ,. _ A." 

■ ■• «•■-. *.ut..r. — i.^ l. j "V jiuim.m,ui visiujiaix^s m mc nicaii lime" 

to failure and standard deviation presented here are applicable only 
when all the n test items are tested until they fail. 
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a) The point estimate of the mean time to failure is: 



fflr 



n s t 



i = l 



b) The point estimate of the standard deviation is 



1=1 



* = IL ( h ~ mw ) _ a In 2 h z - ( £ h ) 2 

a A/ i^l -\/ i=l i = l 

V n - 1 V n(n - 1 ) 



7.2.1.2 Confidence intervals — The confidence limits for the true 
mean time to failure, wf, with a confidence level of 90 percent are 
shown below [ see Table 4 for values of the Student's distribution 

t* ( v )]. 



TABLE 4 STUDENT / DISTRIBUTION 



Degrees of 
Freedom v 



(v) 



'0-95 



Degrees of 
Freedom v 



'0-90 
(") 



^0-95 

(v) 



2 


1-89 


2-92 


14 


1-34 


3 


164 


2-35 


16 


1-34 


4 


1-53 


213 


18 


1-33 


5 


1-48 


2-02 


20 


1-33 


6 


F44 


1-94 


25 


1-32 


7 


1 41 


1-89 


30 


1-31 


8 


1-40 


1-86 


40 


1-30 


9 


1-38 


1-83 


60 


1-30 


10 


1-37 


1 81 


100 


1-29 


12 


1-36 


1.78 


oc 


1'28 



1-76 



75 
73 
72 
71 
70 
68 
67 
65 



1,64 



Note — Linear interpolation for intermediate v values is sufficiently accurate. 
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a) One-sided confidence interval, lower limit: 

a 



777F > m F — t . (j0 (n — 1 ) 

b) Two-sided confidence interval: 

a 



An 



™? — t . 95 ( n — 1 ) ° < m F < wf + t . 9 - o ("—!) — 



An 



'7.2.2 Graphical Method 



















I 




















99- 




































§«; - 


















G£ . 


















&£ - 
















/' 


















70- 




























/^l 


> 


, 


5 0- 
















' 
















I 


"> n- 
















1 










&S 








1 




















10- 
5- 


































2 - 


















1 - 




































































i 


















. . 1 




R 


JO 1200 K00 ' 
ELEVANT TFST TIMF FOR F&HH TIMF 


I 1 

A 
„ ex =195 


1800 

















In example, w=10, r=& 

Failures occurred times in hours as plotted m?— 1 500 

Fig. 5 Example of use of Normal Probability Paper 
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7.2.2.0 The graphical method is based on the normal probability- 
paper. This paper is so designed that the cumulative distribution function 
of times to failure, which follow a normal distribution, will become a. 
straight line, see Fig. 5. It is applicable also in the case that not all 
items have failed. 

7.2.2.1 The observed h values are plotted at the corresponding 

median rank values, /> 5 o ( t\ ), on the F( t ) scale. Median rank values 

t — 0'3 
may be found in Table 1 or approximated by the ratio — Z~~n T 4 - wnere 

n is the number of test items. The example of Fig- 5 assumes a test 
sample w = 10 with the test, terminated after r = 8 failures. In such 
a case, where the test is terminated before all items have failed, only the 
median ranks up to and including i — r are plotted- 

In drawing the line the points between the probabalities 20 percent 
and 80 percent should principally determine the line. 

7.2.2.2 The estimate of the mean time to failure, bif, is equal to- 
the t value of the point on the straight line, which has an ordinate 
equal to 50percent. 

7.2.2.3 The estimate of the standard deviation, cr, is equal to the 
difference between the t value of the point on the straight line which 
has an ordinate equal to 84 percent and the estimated mean time to> 
failure. 

7.2.2.4 The fitted straight line may also be used to estimate the 
fraction of the population below or above a given time to failure. 

8. SUCCESS RATIO 

8.0 The success ratio is the probability that an item will perform a 
required function or that a trial will be successful under stated conditions. 

When testing for success ratio, the test items or trials are to be 
classified as either failed or non-fai!ed. Reusable devices may be repaired' 
between successive trails, provided that the state and performance are 
the same at the start of all trials. The point estimate ( observed success 
ratio ) and the confidence interval for the success ration are independent 
of how the test is terminated at a fixed number of trials, a fixed 
number of successes or a fixed number of failures. 

8.1 Point Estimate — The point estimate of success ratio is equal to the- 
number of successes divided by the total number of trials or test items: 

n — r 

e = 

n 

8.2 Confidence Intervals 

8.2.0 The confidence limits for the true success ratio R with a 
confidence level of 90 percent, based on the binomial distribution,- are 
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shown below. (See Tables 5 and 6 for values of the F distribution for 
0*9 and 95 fractiles. ) 

TABLE 5 F DISTRIBUTION — 9 FRACTILES 



v s \ v i ' 2 



JO 



20 



30 



40 



60 



120 



2 
4 
6 
8 

10 
12 
14 
16 
18 
20 
30 
40 
60 
120 



9-00 

4-32 
3-46 
311 
292 
281 
2-73 
2-67 
262 
2-59 
2-49 
2 44 
239 
2-35 
230 



^•24 
4 11 
3-18 
2 81 
2-61 
2-48 
239 
233 
229 
2 25 
214 
209 
204 
1 99 
1-94 



9-33 

401 
3-05 
2 67 
2-46 
2-33 
2 24 
21 8 
2 13 
209 
1-98 
1-93 
1-87 
1-82 
1-77 



9-37 
3-95 
298 
259 
2-38 
2-24 
215 
209 
2 4 
2 00 
1-88 
1-83 
1-77 
1-72 
1-67 



9 39 
3-92 
2-94 
2-54 
2-32 
2 19 
2 10 
203 
1-98 
1 94 
1-82 
1-76 
1-71 
1-65 
1-60 



944 



■84 
•84 
■42 
■20 
06 
96 
•89 
•84 
•79 
•67 
■61 
•54 
•48 
1-42 



946 


9-47 


9-47 


9-48 


9-49 


3-82 


3-80 


3-79 


3 78 


3 96 


2-80 


2-78 


2-76 


2-74 


2-72 


2-38 


236 


2-34 


2-32 


2 29 


2-16 


2-13 


211 


208 


206 


201 


1-99 


1-96 


1 93 


1 90 


1 91 


1-89 


1-86 


1-83 


1-80 


1-84 


1-81 


1 78 


1-75 


1-72 


1-78 


1-75 


1-72 


1-69 


1-66 


1-74 


1 71 


1-68 


1 64 


1 61 


1-61 


1-57 


1-54 


1-50 


1-46 


1-54 


1-51 


1 47 


1-42 


1-38 


1 -48 


1-44 


140 


1-35 


1-29 


1-41 


1-37 


1 -32 


1-26 


119 


1-34 


1-30 


1 !-2* 


117 


1-00 



Note — 


■ Linear interpolation for intermediate m values is sufficiently accurate. 






TABLE 6 F DISTRIBUTION — 95 FRACTILES 


r.W, 


2 4 6 8 10 20 30 40 60 120 


CO 



2 

4 

6 

8 

10 

12 

14 

16 

18 

20 

30 

40 

60 

120 



1900 
6-94 
514 
4-46 
4-10 
3-89 
3-74 
363 
3-55 
3-49 
3-32 
3-23 
3 15 
3-07 
300 



1920 
6 39 

4-53 
3-84 
348 
3-26 
311 
3 01 
293 
2-87 
269 
2-61 
2-53 
245 
237 



1930 
616 

4-28 
3-58 
322 
3-00 
2 85 
274 
266 
260 
2-42 
2-34 
225 
218 
2-10 



40 
04 
15 
44 
07 
85 
70 
59 
51 
45 
27 
18 
10 
02 
94 



1940 
5-96 

406 
3 35 
2-98 
2-75 
2-60 
249 
2-41 
2-35 
2-16 
2 08 
1-99 
1-91 
1-83 



19-40 
5-80 
3 87 
315 
2-77 
254 
2-39 
228 
219 
2-12 
1-93 
1-84 
1-75 
1-66 
1-57 



1950 
5-75 
3-81 
3 08 
270 
2-47 
2-31 
2-19 
211 
2-04 
1-84 
1-74 
1-65 
1-41 
1-46 



19-50 
572 
3-77 
304 
266 
243 
227 
215 
206 
1-99 
1-79 
1-69 
1-59 
1-37 
1-39 



19-50 
569 
3-74 
3-01 
262 
238 
2-22 
2-11 
2-02 
1-95 
J -74 
1-64 
1-53 
1-32 
1-32 



1950 
5-66 
3-70 
297 
2-58 
2-34 
2 18 
2.06 
I 97 
1-90 
1 68 
1-58 
1-47 
1-26 
1-22 



•50 
•63 
•67 
93 
■54- 
■30 
13 
•01 
•92 
•84 
•62 
•51 
•39 
•19 
■00 



Note — Linear interpolation for intermediate v values is sufficiently accurate. 
8.2.1 One-sided confidence interval, lower limit: 



R > 

with 



n — r 



n — r + ( r + 1) F B(i ( r u v 2 ) 
vi = 2(r + l) and v z = 2-(n 
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8.2.2 Two-sided confidence interals: 
n — r 



<R< 



( n — r + 1 ). F n ( v u v 8 ) 



n - r + ( r -\- I ) F n ( v„ v a ) ^ AX ^ r +"( /*- r + 1 ) F 95 ( v u v % ) 

with vj = 2 ( r + 1 ) and v. 2 = 2 ( n — > )at lower limits, and 

^ = 2(rc— r + 1) and v- 2 = 2r at upper limits 

F.3 Use of Nomographs 

8.3.1 Nomographs are presented for determining one-sided and two- 
sided confidence intervals of success ratio, see Fig. 6 and 7. The 
nomographs are based on the total number of test items or trials, n, and 
the total number of observed failures, r. 



o 
o 




0-0001 0-00) 0-01 0-1 

FAILURE RATIO, r/n ^> 

Fig. 6 
8.3*2 Enter the a^Wicable nomo^ra^h on the abscissa at -the rhi value. 
Go vertically up to the curves labelled with the observed number of 
failures, r. Then go horizontally to the ordinate and to read the relevant 
failure ratio limits, Q. The confidence limits of the success ratio are the 
complements, that is, R = 1 — Q. 



24 



IS : 8161 ( Part 4 )- 1985 




0-0001 
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WiiMBER 



0-001 0-01 

FAILURE RATIO, r/n 
FIG. 7 
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h CALENDER 
TIME 



Fig. 8 
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i u nu \t n t v a 

( Clauses 6.1.0 and 6.2 ) 
DETERMINATION OF ACCUMULATED RELEVANT TEST TIME 

A-l. The relevant test time of each individual test item can be measured 
by elapsed time indicators on the test item. In this case the accumulated 
relevant test time Ik at the Ath failure is the sum of indicator readings. 

n 

Tk = £ tk, m 
m = 1 

where 

n = the total number of test items, and 

h m — the indicated relevant test item of number m up to the Ath 
failure in the total number of test items. 

A-2. The accumulated relevant test time T* at a decision point not 
coinciding with a failure is: 

n 

i - = a r m - 
m = 1 

where 

t^ = the indicated test time of item number m up to the decision 
point. 

A-3. If the relevant test time is recorded by other means, the accumu- 
lated relevant test time 7V, at the Ath failure should be calculated by the 
following iterative formula. This includes the accumulated relevant 
test time up to the failure k — 1 and the relevant test times elapsed in 
the interval between the failures k — 1 and k. 

n 
?k = Tli ~ 1 + £ £j, tm, i 
m — 1 

where 

n = the total number of test items, and 

# i . — +k*i T+V-i «fit*Ir\*-1 s\P -hi*1 n-wo rt+ +^o+ ft m a r*f t^of i + pnn mi mKf»r w* 

4- JJL 5 J iliv j til i/viiwu v/x iwiwyuin vvjv univ va ivji itvin jiuuiltvi. i »-* 

after the failure A: — 1 in the total number of test items. 

A-3. The intermissions may be caused by failure k — 1 or by any other 
technical or administrative reasons. The number, j, of intermissions 
may vary from item to item. 
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A-4. Figure 8 explains the numbering of periods of relevant test time. 

A-5. The accumulated relevant test time T* at a decision point not 
coinciding with a failure is in this case: 



T* = Tr + 2 2, t 



m — I 



m, j 



where 



Ti = the accumulated relevant test time up to the latest failure 
before the decision point, and 

t = the jth period of relevant test time item number m after 

the latest failure in the total number of test items. 

A-6. The above formulae are applicable also in case the test items are 
non-repaired equipments, except that relevant test time does not exist 
rafter the first failure of each test item. 
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INDIAN STANDARDS 
ON 

RELIABILITY OF ELECTRONIC AND ELECTRICAL COMPONENTS 
AND EQUIPMENT 

IS: 

1885 ( Part 39 )- 1979 Electrotechnical vocabulary: Part 39 Reliability 
of electronic and electrical items (first revision ) 

2612-1965 Recommendation for type approval and sampling procedures 
for electronic components 

7354 Guide on reliability of electronic and electrical items: 

7354 ( Part 1 )-1975 Part 1 Preliminary reliability considerations 

7354 ( Part 2 )-1975 Part 2 Managerial aspects of reliability 

7354 ( Part 3 )- 1975 Part 3 Presentation of reliability data on 
electronic and electrical components ( or parts ) 

7354 ( Part 4 )- 1974 Part 4 Collection of reliability, availability 
and maintainability data from field performance 

7354 ( Part 5 )- 1975 Part 5 Inclusion of lot-by-lot and periodic 
inspection procedures in specifications for electronic and 
electrical components ( or parts ) 

7354 ( Part 6 )-1983 Part 6 Inclusion of reliability clauses into 
specifications for components ( or parts ) ( first revision ) 

7690-1975 Mathematical guide to the terms and definitions for reliability 
of electronic equipment and components ( or parts ) used therein 

8161 Guide for equipment reliability testings: 

8161 ( Part 1 )-1976 Part 1 Principles and procedures 

8161 ( Part 5 )-1981 Part 5 Compliance test plans for success ratio 

8161 ( Part 6 )-l983 Part 6 Tests for validity of a constant failure 
rate assumption 

8161 ( Part 7 )-1977 Part 7 Compliance test plans for failure rate 
and mean time between failures assuming constant failure rate 

8161 ( Part 11 )-l983 Part 11 Flow chart describing preparations 
for and execution of reliability tests 

9185 Endurance ( life ) test for electronic and electrical components: 
9185 ( Part 1 )-1979 Part 1 Thermal endurance 
9185 ( Part 2 )-1979 Part 2 Mechanical endurance 

9186-1979 Guide for screening of electronic and electrical items 

9692 Guide on maintainability of equipment: 



9692 ( Part 1 )-1980 Part 1 Introduction to maintainability 

9692 ( Part 2 )-1980 Part 2 Maintainability requirements in speci- 
fications and contracts 

9692 ( Part 3 )-198l Part 3 Maintainability programme 

9692 ( Part 5 )-Part 5 Maintainability studies during the designphase 

9692 ( Part 6 )-1983 Part 6 Maintainability verification 

9692 ( Part 7 )-i984 Part 7 Collection, analysis and presentation 
of data related to maintainability 

10139-1982 Presentation of reliability, maintainability and availability 
predictions 

10673-1983 Sampling plans and procedures for inspection by attri- 
butes for electronic items 

11137 ( Part 2 )- 1984 Analysis techniques for system reliability: Part 2 
Procedure for failure mode and effects analysis ( EMEA ) and 
failure mode, effects and criticality analysis ( EMECA ) 



